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Abstract 

We consider a collapsing sphere and discuss its evolution under 
the vanishing expansion scalar in the framework of f(R) gravity. The 
fluid is assumed to be locally anisotropic which evolves adiabatically. 
To study the dynamics of the collapsing fluid, Newtonian and post 
Newtonian regimes are taken into account. The field equations are 
investigated for a well-known f(R) model of the form R + 5R 2 admit- 
ting Schwarzschild solution. The perturbation scheme is used on the 
dynamical equations to explore the instability conditions of expansion- 
free fluid evolution. We conclude that instability conditions depend 
upon pressure anisotropy, energy density and some constraints arising 
from this theory. 
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1 Introduction 

To accommodate observational data with theoretical predictions, one needs 
to introduce the dark energy (DE) contributions in the most successful grav- 
itational theory of General Relativity (GR). On the other hand, modified 
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theories of gravity may provide a cosmological accelerating mechanism with- 
out introducing any extra DE contribution. Modifications of GR by f(R) 
models are able to mimic the standard ACDM cosmological evolution. These 
models may have vacuum solutions with null scalar curvature that allow to 
recover some GR solutions. These may also lead to the existence of some 
new solutions particularly in spherically symmetric scenario. 

In the last decade, some f(R) models were considered to modify GR at 
small scales to explain inflation, e.g., f(R) oc R 2 but failed to explain the late 
time acceleration. A model like f(R) oc ^ was proposed to explain this accel- 
eration but attained no interest due to conflict with solar system tests (Chiba 
2003; Dolgov and Kawasaki 2003). A cosmological viable model needs to 
satisfy the evolution of big-bang nucleosynthesis, radiation and matter dom- 
inated eras. Also, they must provide cosmological perturbations compatible 
with cosmological constraints from cosmic microwave background and large 
scale structure. The problem of cosmological perturbations in this modi- 
fied theory and their consequences have widely been discussed in literature, 
e.g., (Hwang and Noh 2006; Carroll et al. 2006; Carloni et al. 2008; Tsu- 
jikawa et al. 2008). Since the lagrangian R + f(R) is analytic at it! = 0, so 
the Schwarzschild and other important GR solutions (without cosmological 
constant) are also the solutions of f(R) gravity (de la Cruz-Dombrize and 
Dobado 2006). 

It is known that f(R) models reveal black hole solutions as in GR, there- 
fore it is quite natural to discuss the question of black hole features and 
dynamics of the gravitational collapse in this modified theory. When we 
take conformal transformation of f(R) action, it is found that Schwarzschild 
solution is the only static spherically symmetric solution for a model of the 
form R + 5R 2 (Whitt 1984). Also, the uniqueness theorems of spherically 
symmetric solutions for general polynomial action in arbitrary dimension 
were proposed by using conformal transformation (Mignemi and Wiltshire 
1992). Work on black hole solutions in f(R) gravity has been carried out 
by many authors, e.g., (Moon et al. 2011; de la Cruz-Dombriz et al. 2009; 
Mazharimousavi et al. arXiv/1105.3659v3) 

It is found that evolution of expansionfree spherically symmetric distri- 
bution is consistent with the existence of a vacuum cavity within the matter 
distribution (Herrera and Santos 2004; Herrera et al. 2008 ). Skripkin model 
is the first example of expansionfree evolution (Skripkin 1960). This model 
corresponds to evolution of spherically symmetric perfect fluid distribution 
with constant energy density. We have explored this model in f(R) gravity 
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(Sharif and Kausar 2011) and effects of f(R) DE on the dynamics of dissipa- 
tive fluid collapse (Sharif and Kausar 2010; Sharif and Kausar 2011). Also, 
solution of the field equations for Bianchi models and spherically symmetric 
dust solution are obtained in this theory (Sharif and Kausar 2011; Sharif and 
Kausar 2011; Sharif and Kausar 2011). 

The physical applications of such models lie at the core of astrophysical 
background where a cavity within the fluid distribution is present. It may 
help to study the formation of voids on the cosmological scales (Liddle and 
Wands 1991). Also, the expansionfree fluid evolution causes a blowup of 
shear scalar which results the appearance of a naked singularity (Joshi et 
al. 2002). It is interesting to note that expansionfree models defines the 
two hypersurfaces, one separating the fluid distribution externally from the 
Schwarzschild vacuum solution while other one is the boundary between the 
central Minkowskian cavity and the fluid. Furthermore, expansionfree models 
require anisotropy in pressure and inhomogeneity in energy density. 

The existence of stellar model can be assured if it is stable against fluctua- 
tions. The problem of dynamical instability is closely related to the structure 
formation and evolution of self-gravitating objects. A pioneer work in this di- 
rection was done by Chandrasekhar (Chandrasekhar 1964). Afterwards, this 
issue has been investigated by many authors for adiabatic, non-adiabatic, 
anisotropic and shearing viscous fluids (Herrera et al. 1989; Chan et al. 
1989; Chan et al. 1994; Chan et al 1993; Herrera et al, arXiv/1010.1518). 
The dynamical instability of collapsing fluids can be well discussed in term 
of adiabatic index T. It measures the variation of pressure with respect to 
a given variation of energy density and defines the range of instability. In 
GR, the dynamical instability is independent of this T factor (Herrera et al, 
arXiv/1010.1518) in the Newtonian and post Newtonian regimes. In a recent 
paper (Sharif and Kausar 2011), we have investigated the dynamical insta- 
bility of expansionfree fluid evolution in f(R) gravity by assuming all the 
higher order curvature terms on the matter side. This provides the effects of 
f(R) DE on the instability conditions. 

In this paper, we treat the field equation in the usual fourth order form 
and discuss instability conditions of expansionfree fluid collapse in metric 
f(R) theory. Here we use Ricci scalar to evaluate time dependent part of the 
perturbed quantities. We consider locally anisotropic fluid distribution inside 
the collapsing sphere. The format of the paper is as follows. In section 2, we 
formulate the field equations and dynamical equations for spherically sym- 
metric spacetime in f(R) gravity by taking locally anisotropy in the pressure. 
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Section 3 is devoted to study the perturbation scheme on physical quantities 
and the metric coefficients by assuming that fluid initially is in static equi- 
librium. A well-known physical f(R) model is considered and perturbed on 
the same pattern. In section 4, Newtonian and post Newtonian regimes are 
taken into account. Discussion is made on the dynamical instability condi- 
tions of the expansionfree fluid evolution in the Newtonian approximation. 
The last section 5 concludes the main results of the paper. 

2 Field Equations and Dynamical Equations 

We take spherically symmetric distribution of anisotropic collapsing fluid in 
co-moving coordinates. The line element is given by 

ds 2 = A 2 (t } r)dt 2 - B 2 (t, r)dr 2 - C 2 (t, r)(d6 2 + sin 2 6d(j) 2 ). (2.1) 

This yields two radii: one is the areal radius C (r, t) measuring the radius from 
spherical surface while other is the proper radius found from J B(t,r)dr. 
Also, there is a hypersurface H* 6 -* separating interior spacetime from the 
exterior one. The energy-momentum tensor for locally anisotropic fluid can 
be written as 

T a p = (p + p±)u a U(3 -p±g al 3 + (p r -p±)XaX(3, (pc,p = 0, 1,2,3), (2.2) 

where p is the energy density, p± the tangential pressure, p r the radial pres- 
sure, u a the four-velocity of the fluid and Xa is the unit four- vector along the 
radial direction. These quantities satisfy the relations 

u a u a = l } x a Xa = -l, X a u a = (2.3) 

which are obtained from the following definitions in co-moving coordinates 

u a = A" 1 ^, X a = B~ l 5«. (2.4) 

The acceleration a a and the expansion G are given by 

a a = Ua-pu 13 , 9 = u? a . (2.5) 

Using Eqs.fl23D and (TO> [t follows that 
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where dot and prime represent derivatives with respect to t and r respectively. 
The Einstein-Hilbert action in GR is 

SeH = 1 LJ d * x ^ R ' ( 2 - 7 ) 

where k is the coupling constant, R is the Ricci scalar and g is the metric 
tensor. For f(R) theory of gravity, it is modified as follows 

Smodif = ^f d 4 x^f(R), (2.8) 

where f(R) is the general function of the Ricci scalar. Varying this action 
with respect to the metric tensor, we obtain the following field equations 

F(R)R a p - ]^f{R)g a p - V a V p F{R) + g aP DF{R) = kT^, (2.9) 

where F(R) = df(R)/dR. For spherically symmetric spacetime, these equa- 
tions become 
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The Ricci scalar curvature is given by 
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For the exterior spacetime to S^ e \ we take the Schwarzschild spacetime in 
the form 



2M 



dv 2 + 2drdv - r 2 {d9 2 + sin 2 9d(f) 2 ), 



(2.15) 



where M represents the total mass and v is the retarded time. 

In order to study the properties of collapsing process, we formulate the 
dynamical equations. For this purpose, we use the Misner-Sharp mass func- 
tion defined as follows (Misner and Sharp 1964) 
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m(t, r) = -(1 + g^C,C v ) = -[! + — -—]. 



(2.16) 



This equation provides the total energy inside a spherical body of radius " C" . 
From the continuity of the first and second fundamental forms, the match- 
ing of the adiabatic sphere to the Schwarzschild spacetime on the boundary 
surface, S^, yields 

M S =m(t,r). (2.17) 
The proper time and radial derivatives are given by 

D =L1 D — — — 
T Adf C C'dr 7 



(2.18) 



where C is the areal radius of the spherical surface. The velocity of the 
collapsing fluid is defined as 



U = D r C 



C 
A 



(2.19) 



which is always negative. Using this expression, Eq. fl2.16p implies that 
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The dynamical equations can be obtained from the non-trivial contracted 
components of the Bianchi identities. Consider the following two equations 

J- -ft U-a — U, l a Xa 



which yield 







B C 

P+(P+Pr)^ + 2{ P + P± )- = 0, 

A' C 

Pr + (p+Pr)^ + 2(p r -p ± )— = 0. 



(2.21) 

(2.22) 
(2.23) 



3 The f(R) Model and Perturbation Scheme 

In this section, we consider a particular f(R) model and apply perturbation 
scheme on all the above equations. Consider the following well-known f(R) 
model 

f(R) = R + 5R 2 , (3.1) 

where 5 is any positive real number. The stability criteria for this model 
corresponds to / (R) > 0. For 5 = 0, GR is recovered in which black holes are 
stable classically but not quantum mechanically due to Hawking radiations. 
Since such features are also found in f(R) gravity, hence the classical stability 
condition for the Schwarzschild black hole can be enunciated as / (R) > 
(Sotiriou and Faraoni 2010). 

The purpose of introducing perturbation scheme is to analyze the insta- 
bility conditions of the dynamical equation. We assume that initially all the 
quantities have only radial dependence, i.e., fluid is in static equilibrium. 
Then these quantities have time dependence as well in their perturbation, 
i.e., 

(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
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Also, the Ricci scalar in f(R) model leads to 



R(t,r) = Ro(r) + eT(t)e(r), 
f(R) = R (l + 25R ) + eT(t)e(r)(l + 25R ) 
F(R) = (l + 25R ) + 2e5T(t)e(r), 



(3.10) 
(3.11) 
(3.12) 



where < e C 1. We choose radial part of the areal radius as the Schwarzschild 
coordinate, i.e., Coir) = r. 

Using the above values in fj2.10|) - fl2. 13j) . the static configuration turns out 
to be 
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The static part of the Ricci scalar takes the form 



Mr) = ~bs 



A" 



A'oB'o 



2B' 



Bl 



(3.16) 



Ao A Q B Q B r r 2 r 
Using the perturbed quantities given in Eqs. fl3.2p - fl3.12p along with f )3.13p - 
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(I3.15P in the field equations, it follows that 
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We can write the perturbed configuration of the Ricci scalar curvature, by 
using Eq. fl3.16p . as follows 
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This equation can also be written as 
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For the sake of instability region, we assume that all the functions involved 
in the above equation are such that a remains positive. Consequently, the 
solution of Eq.f l3.22p becomes 



T(t) 



(3.24) 



Here we assume that the system starts collapsing at t — —00 such that 
T(— 00) = 0, keeping it in static position. Afterwards, it goes on collapsing 
with the increase of t. 

Applying static and non-static perturbation scheme respectively on dy- 
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namical equations ( I2.22p and (12.23 p . we have 
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In order to relate p and p r for the static spherically symmetric configuration, 
we may assume an equation of state of Harrison- Wheeler type as follows 
(Chan et al. 239; Wheeler et al. 1965) 
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where V is the adiabatic index. It measures the variation of pressure for a 
given variation of density. We take it constant throughout the fluid evolution. 



4 Expansionfree Newtonian and Post Newto- 
nian Regimes 

Here we assume expansionfree (G = 0) evolution of anisotropic fluid and 
develop dynamical equations. Since the expansion scalar describes the rate 
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of change of small volumes of the fluid, so the expansionfree fluid evolution 
of spherically symmetric distribution is consistent with the formation of a 
vacuum cavity within the fluid. The Minkowski spacetime is supposed to 
present inside the cavity. Using Eq. fl3.9p in (|2.6p . we have 

- f ( b 2c\ j . 

© = t- "5- + - ■ 



A \B r 
The expansionfree condition (G = 0) implies that 

4- = -2^. (4.2) 

Before implementing the above result on the last section, let us identify the 
terms differentiating Newtonian (N), post Newtonian (pN) and post post 
Newtonian (ppN) regimes. These terms will be considered to develop the 
dynamical equations which help to understand the instability conditions of 
the expansionfree fluid evolution. For N approximation, we assume 

Po > Pro, Po > Rlo- (4.3) 
For the metric coefficients expanded up to pN approximation, we take 

Ac=l-^, B = lA (4.4) 

where G is the gravitational constant and c is the speed of light. Adding 
Eqs. fl3TT3|) and (jSHD and using Eq. fl3T29l) . it follows that 

-Aq _ Kr 3 ( P o + p r0 )(l + 25Rp) + 25R' 'r 2 (r - 2m ) + 2m (l + 26R + 6R' r) 
A ~ 2r(r-2m )(l + 2SR + SR' r) 

(4.5) 

Inserting this equation in (13.25j) . we have first dynamical equation in rela- 
tivistic units as follows 
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In view of dimensional analysis, this equation can be written in c.g.s. units 
as follows 

Prt = 2rc-\r - 2G<r*mo)(l + 2^ + Mfr) [GkT (P ° + C ^o)(l + 2^o) 
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(4.8) 

It is worth mentioning here that the order of c differentiates the terms be- 
longing to N, pN and ppN regimes, i.e., 

terms of order c° correspond to N-approximation, (4-9) 
terms of order c~ 2 correspond to pN-approximation, (4-10) 
terms of order c -4 correspond to ppN-approximation. (4.11) 

Simplification of Eq. (14. 8p yields the following table 

Table 1. N, pN and ppN Terms 
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mlp r0 5R ; , p Q ml5R' Q \ m p r0 5R'Q, p r om (l + 25Rq + 5R Q ) 
m 2 p r0 (l + 25R + 5R' ), m 3 p (l + 25R + 5R' ) 
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Using the expansionfree condition (14. 2p . Eq.f l3.28p becomes 

p = 2(p r0 -p ±0 ) — . 

r 



(4.12) 



This shows that the perturbed energy density depends on the static pressure 
anisotropy, hence supporting the expansionfree condition. Inserting Eq. (I4.12p 
in (I3.3ip . we have 
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Notice that the expressions of p r and p^ in Eqs. f l4.12]) -f l4.14j) are of ppN 
order approximation. In order to discuss the instability conditions up to pN 
order, we neglect these quantities. 

Substituting the value of p± from Eq.f l3.20p in dynamical equation (I3.27P 
and using Eqs. (13.241) and (14.21) . it follows that 
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(4.16) 
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This shows that the general dependence of radial function affects dynamics 
of collapsing fluid. For the sake of simplicity, we assume 



a(r) = a 4 air, c(r) = c 4 Cir, e(r) = e 4 eir, 



(4.17) 



where the quantities with subscript "0" and "1" are arbitrary constants. 
Using these values in Eq.f l4.16p . it follows that 
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Cp 4 c x r 2(e + e x r)5 



1 + 25Rq 



Inserting Eqs.flSJ, ( 14.141) and ( 14.151) up to pN order (with c = G = 1) in 
the above equation, we obtain 
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(4.19) 



In view of Table 1 and neglecting terms with p r o/po being of ppN order, 
Eq. (l4.19p reduces to the following equation in the Newtonian regime 

n(r + m )po + ftp r o(3air + 2a — ) + 2K,pj_ (a + a i r + — — 8ci) 
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(4.20) 



In general, the instability range depends upon the index T as it measures 
the compressibility of the fluid. However, the above equation is independent 
of r which shows that instability region totally depends upon the pressure 
anisotropy, energy density, chosen f(R) model and arbitrary constants. No- 
tice that independence of V factor indicates that under expansionfree con- 
dition, fluid evolves without being compressed. In this way, the given f(R) 
model shows the consistency of the physical results with expansionfree con- 
dition. 
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In order to satisfy the instability conditions of expansionfree fluids, we 
need to keep all the terms positive in Eq. fl4.20p . Here we assume that all the 
arbitrary constants and dynamical quantities are positive whereas p' rQ < 
shows that pressure decreases during collapsing process. To keep all the 
terms in dynamical equations positive in the Newtonian regime, we need to 
satisfy the following constraints arising from different terms 
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< ai < 



a + a x r + 
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3c 
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> 8ci, 



2a 



3a x r > — , 
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10c 8ci , 
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(4.21) 
(4.22) 
(4.23) 
(4.24) 



Thus the system would be unstable in N-approximation as long as the above 
inequalities are satisfied. The dynamical equation ( 14.19P for pN regime be- 



comes 
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(4.25) 



r 1 + 25R 

It is remarked that only relativistic effects are taken into account at pN 
regime, however, the dependence of instability condition remains the same. 
Also, the index T does not involve, so instability conditions depend on the 
same parameters and constants as that in the N-approximation. 
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5 Summary 



This paper is devoted to investigate the dynamical instability conditions 
of the expansionfree fluid evolution for a particular model (f(R) = R + 
5R 2 ) in f(R) gravity. To make consistency with the physical application 
of the vanishing G, we have considered the locally anisotropic fluid with 
inhomogeneous energy density. Perturbation analysis is used in N and pN 
approximations and terms differentiated in both the regimes depending upon 
the order of c. 

Expansionfree models may help to study the formation of voids. Voids 
are the spongelike structures and occupying 40%-50% volume of the uni- 
verse. Observations suggest very different sizes of the voids, i.e., mini- voids 
(Tikhonov et al. 2006) to super- voids (Rudnick et al. 2007). As concerned 
to the shapes of the voids, they are neither empty nor spherical. However, 
for the sake of investigations, they are usually described as vacuum spherical 
cavities surrounding by the fluid. The assumption of spherically symmetric 
spacetime outside the cavity is justified for cavities with sizes of the order 
of 20 Mpc or smaller, as the observed universe would be inhomogeneous on 
scale less than 150-300 Mpc. 

The adiabatic index F measures the variation of pressure, its value de- 
fines the range of instability. For example, for a Newtonian perfect fluid, 
the system is unstable for T < 4/3. We have found that like GR (Herrera 
et al. arXiv/1010.1518), our results are also independent of the V factor. 
This shows the consistency of expansionfree condition with f(R) gravity as 
this requires that fluid should evolve without compressibility. Moreover, the 
instability range depends upon the radial anisotropy of pressure, energy den- 
sity and some constraints on the constants arising from the positivity of the 
dynamical equation. Equation (14.201) holds the instability requirement as 
long as the inequalities (I4.2ip - fl4.24p are satisfied. 

It is worthwhile to mention here that the chosen f(R) model is the only 
model admitting Schwarzschild solution. We have used the perturbed Ricci 
scalar to get the time dependent part of the perturbed metric coefficients. In 
this way, dynamics of the gravitational collapse comprises the scalar curva- 
ture in its evolution. 
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